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Abstract Pour alternative functions are used for fitting tree height and diameter growth models for mongolian oak. 
(Ouercus mongolica Pisch. et l ure/,.). The data set includes 1250 random trees and 755 dominant trees corning from 510 
temporary plots. The resultsshow that the Richards function is the best model for predicting height, diameter at breast height 
(DBI1) and dominant height front age. fhe average growth curve of dominant height is used as a guide curve for the con¬ 
struction of a site index table which is partially validated using an independent data set. fhe Milscherlieh function is the best 
model for estimating height and dominant height from DBII. 
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Introduction 

Mongolian Q'dk(Qucrcus mongolica Fisch. et Turcz.) 
stands constitute the largest type of the secondary forest 
in Heilongjiang Province. China(Hu. 1979). The 1991 
statistics of the Administrative Bureau of Forest Re¬ 
sources for Heilongjiang Province indicated an esti¬ 
mated area of 514977 hnr for mongolian oak stands 
and a volume of 32618000 m’ for this species. The 
study of growth and yield of this species is important in 
order to elaborate management guidelines, but such 
studies are still lacking. 

fhe aim of this paper is to develop average height 
and diameter growth models for mongolian oak. We 
first make a short review of some common theoretical 
growth functions and we recall their mathematical 
properties. We then present the data set and compare 
the results obtained when fitting the different growth 
functions. Finally we construct and test a site index 
table for mongolian oak in Heilongjiang Province. 

Data Collection 

The area under study is the Mudanjiang region in 
Heilongjiang Province. It is located between longitude 
128 51 and 134 5 E and between latitude. 43 40 and 
49 25 N. Elevations range from 400 to 800 m. The 
average annual temperature varies between 2.5 C and 
4 (\ The frost-free period is between 110 and 140 days 


and the mean annual precipitation between 460 and 690 
mm -vr' 1 . 

Data were collected from 510 temporary plots lo¬ 
cated in 10 forestry bureaus of the Mudanjiang Forest 
Administrative Bureau in Heilongjiang Province. The 
plots were sampled according to a systematic design. 
All the plots come from oak dominated secondary high 
forest stands that are approximately even-aged. 

2077 trees were sampled in these plots. 827 of these 
trees were selected among dominant trees whereas the 
1250 others were selected randomly (some but few of 
them were dominant ). 755 of the 827 dominant trees 
were used to tit the models and 72 were used to test 
them. Summary statistics are presented in Table 1. The 
measurements that were carried out for each tree are 
total age(A). diameter at breast height(DBH)and total 
height (H). 


Table I. General statistics for stand attributes of the 510 
sample plots 


Stand variables 

Min. 

Max 

Mean 

Standard 

deviation 

Coefficient 

of variation 

stand age(yr ) 

15 

135 

47.0 

20 96 

0 443 

meaji I)till(cm) 

6 

85 

13.6 

7 0.31 

0.516 

mean heighl(m) 

4 

24 

10.2 

3 127 

0507 

stand volume(m7hm ) 

5 

335 

72.0 

4X.56 

0.667 


Stem analysis were carried out on 10 dominant trees 
that were selected and felled in different plots-corre¬ 
sponding to different age and site classes. These data 
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were intended to validate qualitatively the models de¬ 
rived from the other data . 

Methods 

In recent years, the Richards function has been exten¬ 
sively used to fit tree and stand growth models (Pienaar 
et al. 1973: Causton & Venus 1981; Ito et al. 1985; 
Houllier et Leban 1991). This function was derived 
from a generalization of the von Bertalanffy growth 
function)Richards 1959). The integrated form of Rich¬ 
ards function is: 

}'= a [I-A e.xp(-£/)] c ( 1 ) 

where T is a tree or stand attribute(e.g. height, biomass): 
i is time or age; a , b, k and c are parameters: 

a is the upper asymptotic value of Y and k is related 
to time scale: a >0. k >0; 

h is related to the origin of the growth curve on the 
time axis while c determines the shape of the growth 
curve and the relative location of its inflection point 
along the Taxis. Possible combinations of b and c are(b 
>0 and c>=l) or (b <0 and e<0). 

The mathematical properties of the Richards growth 
function can be obtained by deriving equation(l) (see 
Fig. 1) as follows: 

Let Z be the current growth rate. Zr the current relative 
growth rate and Zm be the total mean increment func¬ 
tion: 


/ = — = ethke o\p(-Av)(l -b e\p(-Ar))‘ 
ill 


Z=kcY((j)' u -1) 

(2) 

/r = kc((j)' u - 1) 

(3) 

Zm = = --(! — hcxp{-kt )) 1 

I t 

(4) 

Let Zmax be the maximum growth rate 
7'„,|, x (resp. Jdenote the time(resp. the T 

which Z-Z nm : 

function, 
value) at 

Z m . K =ak( -)‘ 1 

c 

(5) 

T m. K = -^g(hc) 

(6) 

v c-l t 

ni;i\ “ ‘d I 

(7) 


that are derived from the general Richards equation) 1) 
and that are frequently used in forestry are the follow- 
ing(Li 1988. 1991): 


Total growth function vs. age 



Current and mean annual increment vs. age 



Current annual increment vs. total growth 



Y 

Fig. 1 Different aspects of the Richards growth function 

Y = y ield (= cumulated annual increments ) 

/ =current annual increment (= growth rate ) 

Zm = mean annual increment 

a - 16: b=l: k=0.04:c=2 

•the simplified Richards function obtained with c>! 


Different values of b and c give rise to different func¬ 
tions! Jiang et al.l990a; Millier 1982). Some functions 
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and b= 1 so that Y < „. =0 

) =«(l-exp(-A7))‘ (8) 

•the Mitscherlich function obtained with c=l and 
h >0 and which has no inflection point: 

) = a(\ - hcxp{-ki)) (9) 

•the logistic function obtained with b = -b < 0 and 
c=-l: 


i + /r'expt -hi) 

•the Gompertz function obtained as the limit when c 
tends to infinite: 

) = a exp (-</ exp( ki ) (II) 

where d is a new parameter derived from u. h and k. 

In this study we used these simplified growth func¬ 
tions feq .(8) to (II)] to fit individual diameter and 
height growth data for mongolian oak . We also used 
these functions for predicting individual height from 
DBH . 

The parameters were estimated using the NUN pro¬ 
cedure of SAS package(SAS 1988) and choosing the 
Marquardt method (Marquardt 1963). The models were 
fitted separately for the 1250 random trees and the 755 
dominant trees . For equation (8) [resp.(9)]. we system¬ 
atically tested whether c (resp.b)was different from 1 
we fixed it to 1 and we fitted the model again . 

Results 

Comparison of the different equations 

The statistical results are presented in Table 2 (resp. 
Table 3) for the 1250 randomly selected trees (resp. 
755 dominant trees). The four theoretical growth func¬ 
tions [eq.(8) to (11)] have approximately similar per¬ 
formances according to the following criteria: the 
pseudo-coefficient of determination(R : )and the root 
mean squared error (RMSH). The simplified Richards 
function [eq. (8)] is the best model for predicting indi¬ 
vidual height. DBH and dominant height from age. The 
Mitscherlich function [eq .(9)] is the best model for 
predicting height and dominant height from DBH . 

Construction of the site index table. 

I he proportional melhod(.liang et al ..1990 b, Li et al .. 


1990) also called the Schumachers or the guide curve 
methodfClutter et al .1983) was used to derive a set of 
height growth curve fitted to the simplified Richards 
function was used as the guide curve: 

Hu = 16.04[ I - exp ( -0.0408.4 ) ]' ,,7( ' (12) 


Table 2. Results of nonlinear adjustment of eq.(8) to (11) for 
the 1250 randomly selected trees. 


Dependent independen Iq nation 

variable ! variable 

R : 

RMSI 

Remark* 

H 

A 

(X) 

0.634 

2.17m 

c 1 

H 

A 

(9) 

0 634 

2 17m 

b 1 

ii 

A 

(Id) 

0.629 

2 24M 


ii 

A 

(ID 

0.626 

2 17 


Best function is 

11=15 OX | 

1 - exp ( 

-0 0261 A 1 

1 


Dim 

A 

(X) 

0 8X3 

2 92 cm 


Dim 

A 

(9) 

0 872 

5 12 cm 

h 1 

Dim 

A 

(10) 

O.XXI 

2 95 cm 


duh 

A 

(II) 

0 X60 

5 1X cm 


Best function is: 

DBII - 145.0 | 1 - e\p(-0.002X)A 

I 1 "*'" 


11 

DBII 

(X) 

0 674 

2.07m 


II 

DBII 

(9) 

0 694 

2 oont 


II 

DBH 

(10) 

0.697 

2.0lm 


II 

DBII 

(11) 

0 670 

2 (Wm 


Best function is 

11 = 19,64 ( 

1-0 9442 

cypl-O04X9 DBII)| 


* In equalums(X|| resp(9)| 

when c< resp h) was 

not significantly diflerent 

from 1. its v alue 

was arbitrariK set to 1 and the model was lilted 

again . 

Table 3. Results of nonlinear adjustment of eq.(8)to (11) for 

the 755 dominant trees. 




Dependent Independc liquation 

IT 

RMSI. 

Remark* 

variable nt i 

. at tabic 





II 

A 

IX) 

0 6X7 

2 (Him 


II 

A 

(9) 

■0 684 

2 06m 


11 

A 

( 10) 

0 6X2 

2 oi m 


II 

A 

(II) 

0.654 

2 lorn 


Best tunction is 

II - 16 04 | 

1- e\p(-0()4()XA)|' 



DBH 

A 

(XI 

0 X79 

3 30 cm 


dbii 

A 

(9) 

0 869 

3.43 cm 

b -1 

Dltll 

A 

(10) 

0 878 

3.32 cm 


DBII 

A 

(II) 

0.857 

3.58 cm 


Best function is 

DBII - 152 

0 11 - e\p( 

4) 0028A )j ! 

s-.i 


ii 

DBII 

(X) 

0 621 

2 2X m 


ii 

DBH 

(9) 

0.633 

2 25 m 

h -1 

ii 

DBH 

(10) 

0 624 

2 2o m 


ii 

dbii 

III) 

0 619 

2 29 m 


Best function is 

11 = 19.20| 1 

- 1.00X4 , 

L*\p(-0.()41 

DBIDI 



* In equations (X) Resp. (9)|. when c(resph) was not significantly diffetent 
from I its value was arbitrarily sel to 1 and the model was lined attain 


50 years was selected as the base age and the w idth 
of site index classes was set to 2 m . The site index 
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table of mongolian oak is presented in Table 4 and the associated growth curves are drawn on Figure 2 . 



Fig. 2 Height growth curves of dominant mongolian oaks in Heilongjiang Province: data and theoretical curves 


Table 4. Site index table for mongolian oak: dominant height 


growth for different site classes. 


Aage 

6 

8 

Site class height at 50 years 

10 12 14 16 

18 

in 

0.91 

1 21 

1.52 

1.82 

2 12 

2 43 

2 73 

20 

2 49 

3.32 

4 16 

4 99 

5 82 

6 65 

7 48 

30 

3 97 

5 29 

6.62 

7 94 

927 

10.59 

1 1 91 

40 

5 14 

6 85 

8 57 

1028 

11 99 

13.71 

15 42 

50 

600 

8 00 

10.00 

12.00 

14(H) 

16,00 

18 00 

60 

6.61 

8,8! 

1 1.01 

13.22 

15 42 

17 62 

19 82 

70 

7.03 

9 37 

1 1.71 

14 06 

16.40 

18 74 

21 09 

xo 

7.32 

9.75 

12.19 

14,63 

17 07 

19 51 

21.95 

90 

7 51 

tool 

12.52 

15.02 

17 52 

20 02 

22 53 

too 

7.64 

10 19 

12.73 

15 28 

17 82 

20.37 

22 92 

1 10 

7.73 

10.30 

12 XX 

15 45 

1 8 03 

20.60 

23.18 

120 

7 78 

10 38 

12.97 

15 57 

18.16 

20 76 

23 35 

130 

7 82 

10.43 

13 04 

15.65 

18 26 

20 86 

23.47 

140 

7 85 

10.47 

13 08 

15 70 

18 32 

20.93 

23.55 

150 

7.87 

10 49 

13 11 

15.73 

18.36 

20.98 

23.60 


Validation 

We first validated qualitatively the model by comparing 
graphically the data coming from the stem analysis and 
the site index curves(Figure 3). 

We then tested the accuracy of the average guide curve 
with 72 dominant trees which had not been used for 
constructing the site index table. We applied the follow¬ 
ing methodfsee Lang et ah. 1989). We estimated domi¬ 
nant height Hu* from age with eq. (12).We then fitted 
the following linear model: 


He - a + h Hu* 

( 13 ) 

and got: a=0. 1435. />=().9954, R= 0.81. n=72. average 
relative error=-0.76%. 

We then calculated the following F-statistic: 

n n 

ncr + 2 a(h - I )£ H ¥ *-(/>- I ) 2 £ //*, * 2 


//* ,-a-h^ll^r 

!= I 

where: Hu = height of dominant trees, Hu*= estimated 
dominant height, n = number of sample trees, a and b 
are the parameters of equationf 13). We got: F=0.6996< 
F 27( , (0.05) = 3.1277. The average height growth curve 
is therefore validated and the site index table is suitable 
for forestry practice . 

Conclusion 

The above results show that the Richards equation is 
suitable for describing average height and diameter 
growth as functions of age .These functions were suc¬ 
cessfully fitted to randomly selected trees and to domi¬ 
nant trees . Applying the proportional method, the av¬ 
erage dominant height growth function was then used 
as a guide curve in order to derive a site index table . 

It must however be recalled that this method is very 
sensitive to the sampling design and. more precisely. 
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that its reliability depends on the balance between age 
and site classes. When the data are unbalanced, this 
method may yield biased results(often, the height 
growth of older stands is underestimated because there 
is a lack of good sites in the older part of the sample). 
The validation that we carried out in this paper does not 
address this problem and id therefore only partial. Solv¬ 
ing this problem would require true dynamic data: stem 
analysis, permanent plots or short-term incre- 
ments(HoulIier, 1990). Unfortunately, only few such 
data were available in the context of this study: we just 
had stem analysis data coming from 10 trees. Figure 2 


indicates that the theoretical height growth curves are 
qualitatively consistent with the observed growth 
curves (we only represented the four oldest trees for 
which stem analysis data were available). Clearly, this 
data set should now be completed in order to perform a 
fully satisfying quantitative validation . 

The relationship between DBH and height is not a 
true growth relationship. Nevertheless we also applied 
and fitted the usual growth functions to study this allo- 
metric static relationship and we found that the 
Mitscherlich function is the best model for predicting 
tree height from DBH . 



Height growth of mongolian oak in Heilongjiang Province: theoretical curves compared to stem analysis data. The 4 

oldest trees are represented. 


Fig. 3 
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